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Abstract 

We theoretically investigate spin decoherence of a single nitrogen-vacancy (NV) center in di- 
Q_| amond. Using the spin coherent state P-representation method, coherence evolution of the NV 

center surrounded by nitrogen electron spins (N) is simulated. We find that spin decoherence time 
as well as free-induction decay of the NV center depend on the spatial configuration of N spins. 
Both the spin decoherence rate (I/T2) and dephasing rate (I/T2 ) of the NV center increase linearly 
with the concentration of the N spins. Using the P-representation method, we also demonstrate 
extracting noise spectrum of the N spin bath, which will provide promising pathways for designing 
an optimum pulse sequence to suppress the decoherence in diamond. 

PACS numbers: 76.30.Mi, 03.67.Pp, 03.65.Yz, 76.30.-v 
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Nitrogen-vacancy (NV) impurity centers in diamond have been investigated extensively 



for fundamental research |2| and potential applications of quantum information process- 
ing devices j3-5|, and a high- precision room-temperature magnetic sensor Is, 7]. Coherent 
properties of NV centers in diamond play a crucial rule in the applications. Spin deco- 
herence is due to couplings to noisy environments. For a NV center in diamond, major 
noise sources are paramagnetic impurities and nuclear spins. Diamond crystals containing 
> 10 ppm of substitutional single-nitrogen (N) impurities are called type-lb diamonds and 
the spin decoherence is caused by fluctuating N spin baths js-lO]. On the other hand, in 
diamond containing much less N spins, the spin decoherence time T2 is much longer and 
the decoherence is often limited by couplings to ^^C nuclear spins (~1.1 % natural abun- 
dance) BQ, 12|. 



Understanding dynamics of an electron spin system and surrounding electron and nuclear 
spin baths has been a long-standing problem for theoretical investigation on electron spin 
resonance (ESR) in solids. The lineshape and width of continuous ESR spectrum as well 
as pulsed ESR signals including free induction decay (FID) and spin echo (SE) decay were 
successfully described by approximating spin baths by stochastic noise fieHs QQ. In order 



to realize spin-based quantum bits (qubits) in solid-state systems 



17 



investigation of 



spin baths is becoming more critical to understand decoherence. Methods like dynamical 



decoupling 



19 



- |2l| has been studied to suppress the decoherence in knowledge of noise 



spectrum of surrounding spin baths 22 



23| . In semiconductor quantum dots 24j . silicon 25| . 



and diamond with very low concentration of N spins [Ul, |12|, |26| , nuclear spins are the major 
source of the decoherence in the system. The dipolar coupling between nuclear bath spins 
is much weaker than the hyperfine (HF) coupling between an electron (central) spin qubit 
and bath nuclear spins (the system-bath coupling). In such cases, microscopic treatments o: 



nuclear spin baths have been thoroughly studied using the cluster expansion method |27l-l31| 
in which the intra-bath coupling is treated as a perturbation to the system-bath coupling. 
On the other hand, for electron spin baths where the strength of the intra-bath coupling is 
comparable to that of the system-bath coupling (the strong intra-bath coupling regime), it is 



often challenging to implement the cluster expansion method [32|, |33| . In type-lb diamond, 
the dipolar coupling strength between the NV center and the N spins is of the same order of 
that between N spins. Although it has been shown that experimentally observed FID and 
SE signals of a single NV center have been well described by a classical stochastic noise, the 
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Ornstein-Uhlenbeck (0-U) process 10|, 124], spin dynamics of a NV center and surrounding 



N spins has not been fully understooc 
Variations of the P-representations 



yet. 
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37l | have been successfully used to formulate many- 
body problems in quantum optics, e.g. squeezing in optical solitons jssl. The spin coher- 
ent state P-representation, as a variation of the time-dependent mean-field method, has 
been proposed to apply to spin-based qubit systems targeting at decoherence of the central 



spin |39|, |40|. The wavefunction of the whole system is represented in the basis of direct 
product of wavefunctions for individual spins. Equations of motion for all spins are spe- 
cially tailored to achieve a close approximation to quantum dynamics of the central spin. 
Simulation of FID and SE for a nuclear spin system has been demonstrated using the spin 



coherent P-representation method 



40| 



In this paper, we theoretically study spin decoherence of a single NV center in diamond 
with the N electron spin bath. We employ the P-representation method to simulate spin 
dynamics of the NV center and surrounding N spins. Simulated SE decays agree with the 
decays due to a classical noise field described by the 0-U process. Our simulation shows 
that the decay rate of both FID (I/T2*) and SE (I/T2) depend linearly on the concentration 
of N spins in the range from 1 to 100 ppm. Noise spectrum of N spins is also extracted using 
the P-representation method. The simulated noise spectrum is in good agreement with the 
noise spectrum of the 0-U process in the range of high frequencies with which the spectrum 
of the SE sequence overlaps significantly. 

We consider a NV center in diamond {S = 1) under a static magnetic field Bq applied 
along the N-V axis (denoted as the 2;— axis). The Hamiltonian of the NV center system is 

Hs = DiS^f + ^oBoS'o (1) 

where D = 2.87 GHz is the zero-field splitting due to the axial crystal field, 70 is the 
gyromagnetic ratio for the NV center, and Sq is the electron spin of the NV center. The 
second term corresponds to the Zeeman energy of Sq. The HF coupling between 5*0 and 
the nitrogen nuclear spin in the NV center is not considered. The degeneracy between 
rriz = +1 and = — 1 states is lifted by the external magnetic field Bq. We consider only 
the rriz = — 1 and = transition to treat the NV center as a two-level system, sq. The 
Hamiltonian for an individual N spin [S = 1/2) is 

H, = ^BqSI + A^Sm (2) 



where 7 is the gyromagnetic ratio for the N spin and the HF couphng constant is Ai = 114 
MHz for N spins delocahzed along the [111] axis and 86 MHz along the [111], [111] and [111] 



axes 



411 • here consider cases where Bq is not equal to 514 G, therefore there is a large 
mismatch in the transition energy of NV and N spins which suppresses the fiip-fiop process 
between NV and N spins greatly. In the rotating frame with the precession frequency of 
the N spins, '-/BoSk, and with the NV center, {D — 7o-Bo)so, the Hamiltonian for the dipolar 



interaction between the NV and the N spins (the system-bath coupling) is given by 10|, 

i/5B = (.g-l/2)5^A,5,^ (3) 

A: 

where Ak = [1 — 3{nl)'^]ak and = h'-fo'-f/rl is the coupling constant. The Hamiltonian for 
the dipolar coupling between N bath spins Sj and Sk is given by 

He = $^c,,,[l - 3{nl,nS]S^, - \{SlS~ + SrS^)] (4) 

where Cj^k = ^l^/f^^^k o^^y secular terms are considered. 

Fig. [1] shows a spatial configuration of a NV center (Sq) and N spin (Sk) bath used in our 
simulation. The NV center spin is located at the center of a cube where the cube consists of 
a tetrahedral diamond unit cell with a = 3.567A of the lattice constant. The number of the 
unit cells in the cube is {Na is the number of unit cell along the x,y and z-axes). N spins 
are randomly distributed on the lattice sites. Na and the number of N spins (typically 80 
~ 100 N spins are used) are adjusted according to the concentration of N spins (/). In the 
simulation, M initial vectors are sampled to represent unpolarized N spin bath corresponding 
to a high temperature limit of the bath, T ^ h^yBo/ks, (Ro™'',Ri"*\ ■ ■ ■ ,Rj™''') and m = 
1, . . . M, are prepared where J is the number of N spins and R/j. = {6k, (pk) is a classical 
vector on the Bloch sphere of spin k (the k=0 spin is a NV center). Time evolution of the 
vector, Rfc(t), is calculated according to a set of equations of motion R^ = x R^ where 
Bfc is a local magnetic field for spin k induced by other spins in the system. Coherence at 
time t is given by calculating (sg(t)) = J2m=i R^'^"^K't) (see the supplemental material for 
details). 

Simulated coherence are shown in Fig. [H^b). The coherence decays as a function of time, 
and the decay corresponds to FID signals of a NV center. As shown in Fig. [^b), simulated 
FIDs are fit well by the fundamental Gaussian function. The Gaussian shape of the FID 
agrees with that of FID signals experimentally observed from a single NV center in type-lb 




1 1.5 2 

Time (f^s) 



FIG. 1: (a) A Schematic of a spatial configuration for a NV and N spins. The cube consists 
of the diamond lattice. A NV center (red large sphere) is located at (0, 0, 0) and N spins (blue 
small spheres) are located randomly on diamond lattice sites, (b) Simulated FID signals of a NV 
center with three different N spin bath configurations. The concentration of N spins is / = 10 ppm 
for all cases. Blue crosses, green squares and red circles are simulated results, and lines are fit to 
exp[— (t/T2 )^]. Obtained T2 = 0.97, 0.69, and 0.61 fis from the fit are in good agreement with 
calculated values directly from the N spin bath configuration where T2 = V^/b are 0.95, 0.67, and 
0.62 fis respectively. The inset shows a histogram of b for the 90 instances with / = 10 ppm. Black 
solid line represents the expected probability distribution of b (see the supplemental material for 
details). 



diamond 



10|. The root-mean-square of the spin-bath couphng b = lyXlj quantifies the 
FID for the NV spin, i.e. Tg* = a/2/6. Thus the FID time depends on the configuration 
of the local bath spins around the NV center. Fig. [I](b) shows the simulated FID signals 
for the single NV center with different N bath configurations for / = 10 ppm of the N 
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FIG. 2: (a) Simulated SE and FID signals with / = 10 ppm. Circles and crosses are simulation 
results, and a red solid line shows a fit to Eq. ([5]). With b = 1.44 fj,s~^ from the FID data, we 
obtained tc = 2.78 /is. (b) Simulated SE signals for 90 instances of the N spin bath configuration 
with / = 10 ppm. Blue circles are simulation results, and green lines are fits to Eq. The inset 
shows a histogram of a for the 90 instances when SE signals are fitted by exp[—{t/T2)°^]. 

spin concentration. We found that the values of b obtained from Tg* agree with b directly 
calculated from the configuration of the N bath spins. Inset of Fig. [U^b) shows a histogram of 
b. The distribution of b agrees with the theoretically expected distribution 42|. To represent 
a typical configuration of the N spin bath, simulated results with any pairs of spins coupled 
much stronger (50 times) than typical coupling strength are excluded. A deviation at high b 
values is due to this exclusion (see the inset of Fig. [Hb)). We obtained Tg* = 0.02 ~ 1.73 fis 
with 0.53 fis of the the mean value for / = 10 ppm of the N spin concentration. 

As shown in Fig. EJj^a), we simulated the time evolution of SE. The rephasing tt— pulse is 
assumed to be perfect and instantaneous in the simulation. SE decay in electron spin baths 
has been described by treating the bath to be a classical noise field where the noise field 
B{t) was modeled by the 0-U process with the correlation function C{t) = {B{0)B{t)) = 



6 



6^ exp(— |t|/rc). tq is the correlation time of the bath, which measures the rate of the flip- 
flop process between the bath spins. The corresponding noise spectrum is Lorentzian with 
power 6^ and the half- width at half-maximum (HWHM) I/tc- SE decay subject to the noise 



due to the 0-U process is given by 



13| 



E{t) = exp[-{bTcY[t/rc - 3 - e"*/^^ + 4e-*/(2-c)]]. (5) 

In the quasi-static limit (brc ^ 1) indicating slow bath dynamics, E{t) = e"*^*^/'-^^^'^) ~ 
exp[— (t/T2)^]. On the other hand, in the motional-narrowing limit {Btc ^ 1), E(t) = 
e-*/fc exp(— t/T2). As shown in Fig. |2]^a), we found good agreement between our simu- 
lation results and Eq. ([5]). Using the value of b determined by the bath configuration and 
confirmed by the FID, we determined tc- Fig. [2]^b) shows 90 simulated SE decays with /=10 
ppm. For many cases, we found brc > 1, and SE decays are well described by exp[— (t/T2)"] 
where the exponent a is typically between 2 ~ 3 as shown in inset of Fig. Mjo). This is in 



consistence with previous experimental results {a ~ 3) [3J, |43|. We also found a few cases 
in the motional-narrowing regime {brc < 1), in which the SE decay shapes are close to a 
single exponential function. 

We examined the concentration dependence of the FID and SE decay in the range from 
1 to 100 ppm. In Fig. [31(a) and (b), I/T2* and I/T2 are shown for the concentrations of 
/ = 1, 5, 10, 50, 100 ppm. For each value of /, 80 ~ 100 random configurations of 
the bath are simulated and the mean of I/T2* and I/T2 with the standard deviation as an 
error bar are shown. We found a linear dependence on / for both I/T2* and I/T2, namely 
1/T; = 0.19/ and I/T2 = 5x 10"^/ corresponding to T2 = 2.03 and 0.2 fis for / = 10 and 100 
ppm respectively. The linear dependence of the concentration has been experimentally seen 
in I/T2 of ensemble N spins in diamond [4^]. In addition, a linear concentration dependence 
of 1 /T2 has been reported in phosphorous donors in silicon 



321]. The linear dependence can 



be understood as stemming from the dipolar nature of both the spin-bath and the intra-bath 
couplings. Consider the case of brc ^ 1. The average spin-bath coupling oc 1/f^ oc /, 
hence 6^ oc where f is the average distance between neighboring spins. Similarly for N 
bath spins, Tc is roughly proportional to average coupling, which is proportional to 1//. 



For brc > 1, the expansion of Eq. [5] to the leading order in t/rc yields I/T2 oc {h^ /tcY^'^ ■, 
therefore I/T2 is proportional to /. 

Finally we simulated noise spectrum of the N spin bath, S{uj). At time t, the local 
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FIG. 3: (a) l/Tg as a function of the N spin concentration. Each cross is the mean value and 
the error bar is the standard deviation of all instances. Black solid line is a fit to a linear function, 
(b) I/T2 as a function of the N spin bath concentration. Each blue circle is the mean value and 
the error bar is the standard deviation of all instances. Blue solid line is a fit to a linear function. 
The red cross and its error bar are the results with the HF interaction between N electron and ^^N 
nuclear spins taken into account. Difference in simulated I/T2 between with and without the HF 
coupling is smaller than the deviation (see the supplemental material for details). 

magnetic field due to N spins at the NV spin is expressed by Bj^yif) ~ Ylj ^o,i(-R|(^)) where 
the sum is taken over all N bath spins. Direct Fourier transform of B]srv{t) renders the noise 
spectrum, S{u). Fig. H] shows the simulated noise spectrum of N spins. Coherence of the NV 
center decays due to couplings to environmental noise of N spins. The amount of the decay 
is determined by overlaps between the noise spectrum and the spectrum of a pulse sequence 
used in ESR measurements. In SE measurement, the power spectrum of the pulse sequence 
is given by | sin^(ci;t/4)/(a;/4)p where t is the total evolution time 23|. As shown in Fig. HI 
with t = T2, significant overlap between the noise and the pulse-sequence spectrum happens 
at high frequencies where the noise spectrum agrees well with the spectrum of the 0-U 
process. This supports the observation of good agreement between our simulated SE signals 
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FIG. 4: Noise spectrum of N spins with / = 10 ppm. Blue circles are the simulated noise spectrum. 
Red solid line shows the noise spectrum of the 0-U process expressed by the Lorentzian function 
with 1/tc = 0.134 /is~^ where tq is extracted from the SE signal. Black dotted line shows the 
power spectrum of SE pulse sequence | sin^ (wt /4) /(a;/4) p with t = T2 = 3.8/is. 

and the analytical solution with the 0-U process (see Eq. [5]). On the other hand, in the 
present case, the dephasing time T2* of the FID signals is determined by e~* -^-^^ ^i'^)'^ where 
T2* does not depend on details of the noise spectrum S{uj). We found excellent agreement 
in between obtained from the noise spectrum and Tg* extracted from simulated FID. 

In summary, we investigated spin decoherence of a single NV center coupling to N elec- 
tron spins using the P-representation method. The P-representation is a modified time- 
dependent mean-field technique that is suitable to simulate coherence evolution of a central 
spin surrounded by spin baths. Simulated results for SE signals are in good agreement with 
analytical expressions based on the 0-U process and with previous experimental results. We 
found that the FID and SE decay rates, l/Tj* and I/T2, depend linearly on the concentra- 
tion of N spins. We also demonstrated simulating the noise spectrum of the electron spin 
bath using the P-representation method. The P-representation method is suitable to study 
decoherence in other spin-based qubit systems even in the strong intra-bath coupling regime, 
e.g. phosphorus donors in isotropically pure silicon (^^Si < 50ppm). This capability will 
facilitate us to find optimum DD sequences to effectively suppress decoherence in spin-based 
qubit systems. 
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SUPPLEMENTAL MATERIAL 
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Spin decoherence and electron spin bath noise of a 
nitrogen-vacancy center in diamond 

Z.-H. Wang and S. Takahashi 



Here we provide details of the method to simulate spin dynamics of a single nitrogen- 
D . vacancy (NV) center and surrounding single-substitutional nitrogen (N) spins, and discussion 



of the hyperfine (HF) coupling between N electron and ^^N nuclear spins and the spin bath 
configuration dependence of the root-mean-square of the spin-bath coupling (6). 



P4 (i) P-representation Method 
-4—^ ■ 

■ Free- induct ion decay (FID) and spin echo (SE) decay signals of the NV center, and noise 

^ . spectrum of N spins in diamond are calculated using the P-representation method. In our 



simulation, a single NV center is positioned at the center of a cube consisting of a diamond 



^ ■ lattice, and N spins are randomly located on diamond lattice sites. The number of the unit 

5^ ! cells in the cube and the number of N spins are adjusted according to the concentration of 

■ N spins (/). We consider the magnetic dipole interaction between all pairs of spins in the 

^ . system to simulate their spin dynamics. In order to represent typical cases, we excluded 

^ ■ bath configurations in which a dipolar coupling strength between any pairs of spins is more 

than 50 times larger than the typical coupling amplitude, iJ^sf / , where g = 2 is the 
^ I electron (7-factor, /i^ is the Bohr magneton, and a = 3.567A is the lattice constant of the 

^ ■ diamond lattice. 

Initially M = 40000 sets of random vectors (r{,"'\ Ri™\ ■ ■ ■ ,^5^^), m = 1,...M 
representing unpolarized initial states p(0) = ^]^\x){x\ ^^=1 Ifc prepared ^ where 
R-fc = (^fci 0fc) represents a classical vector of spin k on the Bloch sphere. Then we compute 
time evolution of the sample vectors according to a set of equations of motion, = x R^, 
where is a local magnetic field for spin k created by all other spins in the system 
For the NV center spin {k = 0), B^ = y/^/'2J2k'jto^k,k'^k'^ where Al^^, is the z com- 
ponent of the dipolar coupling to spin k. The x and y components are ignored because 
of a negligible flip-flop rate between NV and N spins. For bath spins {k = 1,...N), 
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FIG. SI: Simulated SE signals with / = 50 ppm of the N spin concentration, (a) and (b) shows 
SE signals with two difference instances. Simulated results without (with) the HE coupling are 
shown by filled red circles (blue triangles). Solid red (broken blue) lines are corresponding fits for 
the case without (with) the HE coupling. 

= V^/'^{Y^k'j^k^k,k'^k' ~ \^k,Q^P,z) where 13 = x,y,z and the second term, 1/274^^0, is 
resulted from the mapping of the NV center spin S* = 1 to the effective spin of Sq = 1/2 (see 
Eq. (3) in the main text). Coherence at time t is computed as (so(t)) = jj J2m=i R^'^"^Kt) 
(ii) Effect of the HF interaction between N electron and ^'^N nuclear spins 
Most of N electron spins interact with ^"^N nuclear spin (1=1, the natural abundance 
of ^^N is 99.6%) via the HF interaction. The HF interaction splits the energy levels of the 
N electron spin. Because of a large energy mismatch between N electron and nuclear spins, 
the fiip-fiop process between the electron and nuclear spin is largely suppressed, leaving the 
effective HF term in Hamiltonian of the system to be AkS^I^. The HF interaction provides 
a static detuning magnetic field for N electron spins. In the presence of the static field Bo 
along [111], the HF coupling is 114 MHz for the N spin delocalized along the [111] and 
86 MHz along the [111], [111] and [111], therefore the N spin state is split into five states 
with 1:3:4:3:1 population ratio. While N electron spins in the same energy state can fiip-fiop 
due to the dipolar interaction, the fiip-fiop rate between electron spins with different states 
are negligible due to the large energy difference. With the five energy levels for N spins with 
the magnetic field Bq along the [111] axis, the number of N spins at the same energy level 
decreases, therefore the average distance between N spins with the same energy becomes 
larger. This may result in a smaller l/r^, and longer T2 in SE signals. Fig. ISTT a) shows the 
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simulated SE decays with and without the effect of the HF couphng with the same bath 
configuration (the number and location of N spins are same). As shown in Fig. ISlT a). the 
SE signals with and without the HF effect show no pronounced difference and both decays 
fit well with Eq. (5) in the main text. On the other hand, in the other configuration, the 
SE signal with the HF coupling shows a longer T2, as shown in Fig. ISlT b). We simulated 
80 instances for statistical analysis of the HF coupling effect. Fig. 3 in the main text shows 
the effect of the HF coupling with the distribution due to the spin bath configuration for 
/ = 50 ppm. The obtained I/T2 is 3.67 ± 2.55 fis~^ and 3.35 ± 2.59 fis^^ with and without 
the HF coupling. Thus, we found that the effect of the HF is much smaller than the deviation 
due to the spin bath configuration. 

(iii) Spin bath configuration dependence of the spin-bath coupling constant b 
The value of b depends on the configuration of N spins in diamond. The distribution of 
the spin-bath coupling constant b is given by 3|, 
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where F is determined by the decay rate of FID in an ensemble measurement. Fig. IS2I shows 
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FIG. S2: Simulated FID for ensemble of NV spins with / = 10 ppm of the N spin concentration 
(red dots). A black line represents a fit to a single exponential function giving T2 = l/F = 0.5 yus. 



a simulated ensemble FID result of the NV centers to determine F. The ensemble FID 
result is obtained by averaging over FID results of a single NV center from 90 difference 
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configurations. 
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